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We present infrared colours (in the 25-500 yum spectral range) and UV to radio continuum spectral energy distributions of a sample of 51 nearby 
galaxies observed with SPIRE on Herschel. The observed sample includes all morphological classes, from quiescent ellipticals to active starbursts. 
Active galaxies have warmer colour temperatures than normal spirals. In ellipticals hosting a radio galaxy, the far-infrared (FIR) emission is 
dominated by the synchrotron nuclear emission. The colour temperature of the cold dust is higher in quiescent E-SOa than in star-forming systems 
probably because of the different nature of their dust heating sources (evolved stellar populations. X-ray, fast electrons) and dust grain properties. 
In contrast to the colour temperature of the warm dust, the /350//500 index sensitive to the cold dust decreases with star formation and increases 
with metallicity, suggesting an overabundance of cold dust or an emissivity parameter /3 <2 in low metallicity, active systems. 
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1. Introduction 

The energetic output of any extragalactic source can be deter- 
mined by constructing its spectral energy distribution (SED). 
The stellar component emits in the ultraviolet (UV) to near- 
infrared (NIR) domain, young and massive stars dominating 
the ultraviolet UMj and old stars the NIR. Dust, produced by 
the aggregation of metals injected into the interstellar medium 
(ISM) by massive stars through stellar winds and supernovae, 
efficiently absorbs the stellar light, in particular that at short 
wavelengths, and re-emits it in the infrared domain (5/jm-lmm). 
At longer wavelengths, the emission of normal galaxies is gen- 
erally dominated by the loss of energy of relativistic electrons 
accelerated in supernovae remnants (Lequeux 1971; Kennicutt 
1983) (synchrotron emission). Reconstructing SEDs is thus of 
fundamental importance for quantifying the relative contribution 
of the different emitting sources to the bolometric emission of 
galaxies and studying the physical relations between the various 
galaxy components (e.g., interstellar radiation field, metallicity, 
dust and gas content, magnetic fields). In particular, the impor- 
tance of the infrared domain explored by Herschel resides in the 
dust that, by means of the absorption and scattering of UV, op- 
tical and NIR photons, modifies the stellar spectra of galaxies. 



SEDs are thus crucial for quantifying dust extinction and recon- 
structing the intrinsic distribution of the different stellar popula- 
tions within galaxies. Furthermore, fitting infrared SEDs is nec- 
essary for measuring the dust properties such as mass, tempera- 
ture, fraction of PAHs, and hardness of the interstellar radiation 
field (ISRF), all crucial ingredients in the study of the physical 
processes within the ISM (e.g., Draine et al. 2007). 
The interpretation of the infrared SEDs of normal galaxies has 
already been the subject of several studies (Dale & Helou 2002; 
Dale et al. 2001; 2005; 2007; Chary & Elbaz 2001) even within 
the Virgo cluster region (Boselli et al. 1998; 2003). These how- 
ever were generally limited in the infrared to /I < 170 /vm, the 
spectral domain covered by ISO or Spitzer. These new Herschel 
data allow us to extend, for the first time for such a large sam- 
ple of normal galaxies, to the sub-mm domain (A < 500 fim) 
where the emission is dominated by the coldest dust compo- 
nent. This unexplored spectral range is crucial for determining 
the total mass of dust and for an accurate determination of the 
total infrared luminosity. This paper presents the first, Herschel- 
based statistical study of the properties of the SED of a sample of 
nearby, normal galaxies spanning a large range of morphological 
type and luminosity. 



* Herschel is an ES A space observatory with science instruments pro- 
vided by Principal Investigator consortia. It is open for proposals for 
observing time from the worldwide astronomical community. 

' The UV emission of early-type galaxies is however due to an old 
stellar population (UV upturn; e.g. O'Connell 1999; Boselli et al. 2005) 



2. The data 

Galaxies analysed in this work were observed during the 
Herschel (Pilbratt et al. 2010) SPIRE (Griflin et al. 2010) science 
demonstration phase as part of the Herschel Reference Survey 
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(HRS), a guaranteed time key project designed to observe with 
SPIRE a volume-limited, K-band-selected, complete sample of 
nearby galaxies (Boselli et al. 2010), and the Herschel Virgo 
Cluster Survey (HeViCS), an open time key project focused on 
covering 60 sq.deg. of the Virgo cluster with PACS and SPIRE 
(Davies et al. 2010). To these, we added two galaxies of the 
Very Neai'by Galaxy Sample, M81 (Bendo et al. 2010) and M82 
(Roussel et al. 2010). The present sample is thus composed of 
51 objects with photometric data in the three SPIRE bands out 
of which 33 are Virgo members, 13 background, 3 isolated and 
2 (M81 and M82) nearby galaxies. 

The Herschel data used in this work were processed using the 
Level 1 procedures described in Pohlen et al. (2010), fluxes be- 
ing multiplied by a factor of 1.02, 1.05, and 0.94 at 250, 350, 
and 500 //m respectively, to take into account the updated flux 
calibrations. Integrated flux densities were extracted using the 
QPHOT task of IRAF. We assume a conservative uncertainty in 
the flux density of < 30%. This includes uncertainties on the ab- 
solute flux calibration (15 %; Swinyard et al. 2010), uncertain- 
ties introduced by the map-making technique and the possible 
contamination of background objects (on and off source), which 
might affect the flux extraction procedures. Independent obser- 
vations of the galaxies NGC 4438 and NGC 4435 (the first one 
being an interacting system with extended morphology, the sec- 
ond one a point-like early-type source) performed during both 
the HRS and the HeViCS surveys, provided consistent results 
to within 10 %. This should thus be considered as an upper 
limit to the uncertainty introduced by map making and flux ex- 
traction. Absolute flux calibration uncertainties, being system- 
atic in the three bands, do not affect the observed trends in the 
SPIRE colour-colour diagrams. The dataset analysed here in- 
cludes SPIRE data and measurements available at other frequen- 
cies, from UV to radio centimetre. Most of these data were taken 
from the GOLDMine database (Gavazzi et al. 2003). 
Despite this sample not being complete in any sense, and be- 
ing dominated by Virgo cluster galaxies for which perturbations 
induced by the cluster environment may lead to systematic dif- 
ferences in the emission properties relative to similar isolated 
objects even in the FIR spectral domain (Boselli & Gavazzi 
2006; Cortese et al. 2010), this is the first sample observed with 
Herschel that is suitable for a statistical analysis since it consists 
of well-known nearby galaxies spanning a wide range of both 
morphological type and luminosity. 

3. Far infrared colours 

A phenomenological, model-independent technique for quanti- 
fying the spectral properties of galaxies is that of determining 
their colours. To do this, we combine SPIRE and IRAS flux 
densities, the latter being sensitive to the warm dust component. 
Figure [1] shows the IR colours of the sample galaxies. 
A first analysis of Fig. [1] indicates that in star-forming galax- 
ies the flux density ratios /60//500, /25//250, or /100//250 
are strongly correlated with the generally used IRAS colour in- 
dex /60//100 (panels a, b and c). The dynamic range covered 
by /60//500, however, is a factor of about 30 larger than that 
covered by the /60//100 flux density ratio, and is thus a much 
clearer tracer of the average temperature of the bulk of the dust 
component. Starburst galaxies, generally defined as those ob- 
jects with /60//100 > 0.5 (Rowan-Robinson & Crawford 1989) 
have /60//500 spanning from ~ 3 to ~ 30. The prototype star- 
burst galaxy in the local universe M82 has a /60//500 of ~ 26, 
significantly larger than NGC 4491, a starburst in the Virgo clus- 
ter, and VIIIZwl82, a background merging system at z = 0.07. 
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Fig. 1. The infrared colours of our targets. Galaxies are coded 
according to their morphological type: magenta empty circles for 
E-SOa, red filled circles for Sa-Sb, green triangles for Sbc- Scd, 
blue squares for Sd, Im, BCD, and Irr galaxies. The black dotted 
line indicates the colour expected from the Dale & Helou (2002) 
empirical SED, the red long-dashed line those from Chary & 
Elbaz (2001), the blue-short dashed, and the green dashed-dotted 
line the colours of the morphology- and luminosity-dependent 
templates of Boselli et al. (2003). 



Early spirals (Sa-Sb, red filled dots, see Table[D have /60//500 
colours generally colder than Sbc-Scd (green triangles) and Sd, 
Im, BCD, and Irr (blue squares). Early-types with a synchrotron- 
dominated IR emission (M87, M84) are well separated in all 
IRAS-SPIRE or SPIRE colour diagrams with respect to the other 
dust-dominated E-SOa. The remaining early-types have colour 
indices indicating that the cold dust has a higher temperature 
than in star-forming systems. As for the interpretation of the FIR 
properties of the early-type galaxies in the SINGS galaxy sample 
(Draine et al. 2007), higher dust-weighted mean starlight intensi- 
ties can explain the high FIR colour temperatures of E/SO galax- 
ies. However, the relative importance of X-ray heating (Wolfire 
et al. 1995), stochastic heating, heating from fast electrons in the 
hot gas, as in supernovae, and the (unknown) size-distribution of 
dust grains in these environments with low-density ISM needs 
further exploration. 
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Table 1. Median colours with 1 a- standard deviation in the 
colour distribution for different morphological classes 
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Type f60/fl00 f250/f350 f350/f500 f25/f250 flOO/f250 



f60/f5{)0 



O 

2 0, 



E-SOa' 0.37±0.06 

Sa-Sb- 0.34±0.12 

Sbc-Scd{).43±0.1{) 

Sd-In- {).54±0.23 



2.60±0.38 
2.50±0.24 
2.29±0.40 
2.54±0.35 



3.13±0.39 
2.86±0.22 
2.70±0.44 
2.97±0.35 



{).{)9±0.0I 
0.08±0.09 
0.08±0.04 
0.24±0.24 



I.82±0.73 
1.10±0.42 
1.21±0.45 
1.84±0.63 



5.69±2.96 
3.14±2.66 
3.66±2.60 
I().30±9.{)7 



O 
CD 



Notes: 

1: excluding the synchrotron-dominated M84 and M87 

2: excluding the starburst NGC 4491 



The empirical SEDs of Dale & Helou (2002), Chary & Elbaz 
(2001), and Boselli et al. (2003), despite possible uncertain- 
ties in the absolute flux calibration (15 %), only qualitatively 
cover the wide range of infrared colours observed in our sample 
even excluding the radio galaxies M87 and M84, and underpre- 
dict the /250//350 ratio for a given /100//250 ratio (panel d). 
Furthermore, these models do not reproduce the coldest colour 
temperatures observed in the SPIRE colour diagram /350//500 
versus /250//350 (panel f). It is also interesting that even the 
most active galaxies such as M82 and NGC 4491, which are ex- 
pected to be dominated by warm dust heated by the dominant 
starburst, host a cold dust component as traced by the 500 /^m 
emission which is underestimated by models (see panel a). 



4. Spectral energy distributions 

Combining integrated flux densities from UV to radio centimetre 
we constructed the observed SED of the target galaxies. Figure[2] 
shows some examples of UV to radio centimetre SEDs of galax- 
ies according to their morphological type. Figure |2] shows that 
in the elliptical galaxies M87 and M84, the sub-mm domain is 
dominated by synchrotron emission. M87 is a powerful radio 
galaxy (Virgo A), where synchrotron dominates down to ~ 10 
yum (Baes et al. 2010). M84 is a moderately active radio galaxy 
with a luminosity at 20 cm of 2 x lO^-' W Hz"'. In spirals, the 
SPIRE data closely follow a modified black body (13-2) of tem- 
perature T^ 20 K (dashed line) (e.g., Bendo et al. 2003; 2010). 
This however can be taken just as a first order approximation 
since quantitative data in relation to dust can be determined only 
after an accurate SED fitting. To identify the heating sources of 
the emitting dust, we can use any tracer of the hardness of the 
interstellar radiation field. Here we adopt the birthrate parameter 
b, defined as the ratio of the present star formation rate (SFR) 
to the SFR averaged along the life fo of the galaxjQ hence b oc 
SFRh)IM,. Following BoseUi et al. (2001, 2009), SFR is pro- 
portional to the extinction-coiTected UV or Ha flux, and M, to 
the NIR flux. Therefore b is tightly related to the hardness of the 
UV radiation field. Figure [3] shows the relationship between the 
two colour indices /60//100 and /350//500 and the birthrate 
parameter, this last determined for late-type galaxies only. 
Figure [3] shows that the colour index /60//100, sensitive to the 
presence of warm dust, increases with b, indicating that galaxies 
with the warmest dust temperature are those at present most ac- 
tive in star formation {b> 1). In contrast, the temperature of the 
cold dust appears to be anticorrelated with b, indicating an ex- 
cess of the cold dust emission in the most active galaxies. These 
trends with b may be non-universal since they might be related 
to the presence of cluster galaxies, which are characterised by a 
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The birthrate parameter is also called the specific star formation rate 



Fig. 3. The relationship between the infrared colours /60//100 
and /350//500 and the birthrate parameter b (a, b) and the 
metallicity (c, d). Symbols are coded as in Fig. [1] The dashed 
lines give the linear best fits to the data: /60//100 - 0.18+0.04 
log b + 0.41+0.02 (7?=0.69; p=99.99%); /350//500 = - 
0.23+0.10 log b + 2.72+0.06 (/;=0.39; p=95.16%); /60//100 
= -0.44+0.18 [12+log(0/H)] + 4.22+1.55 (7^=0.55; p=96.59%); 
/350//500 = 1.41+0.55 [12+log(0/H)] - 9.56+4.70 (/?=0.55; 
p=97.40%), where R is the correlation coefficient and p is the 
Spearman's probability that the two variables are correlated. 
Best fits with metallicity were performed excluding the 2 out- 
liers M82 and the perturbed system NGC3448 (Arp 205) since 
their uncertain metallicity is probably not representative of nor- 
mal galaxies. 



reduced star formation activity {b ~Q.\) because of their interac- 
tion with the cluster environment (e.g., Boselli & Gavazzi 2006). 
An opposite behaviour of the /60//100 and /350//500 colour 
indices is also present when plotted versus the gas metallicity in- 
dex 12+log(0/H) (determined using the prescriptions of Kewley 
& Ellison (2008) based on the Pettini & Pagel (2004) calibration 
and using mainly the Gavazzi et al. 2004 spectroscopic data), 
i.e., while /60//100 decreases with metallicity (with the excep- 
tion of the starburst M82), /350//500 seems to increase with 
12+log(0/H), with a possible exception for the interacting sys- 
tem NGC 3448 (Arp 205). A similar radial trend with metal- 
licity is also observed for both M99 and Ml 00 (Pohlen et al. 
2010). This result is consistent with the presence of emission at 
A > 850 yum that could originate in < 10 K dust (Galliano et al. 
2005; Galametz et al. 2010; O'Halloran et al. 2010) or dust with 
y6 < 2 (e.g. Bendo et al. 2006), which may be more prominent in 
low metallicity galaxies. A value of /3 < 2 implies an enhanced 
contribution from carbonaceous dust because amorphous hydro- 
carbons have values of/3 in the range 1-1.5. Before attempting to 
determine the origin of this cold dust component, this interesting 
result should be confirmed on a more robust statistical basis. 



5. Conclusions 

Our analysis has enabled us to reach the following conclusions: 
a) the infrared colour index /60//500 is more capable of detect- 
ing a starburst than /60//100. b) Normal galaxies show a grad- 
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ual increase in their dust temperature along the Hubble sequence, 
from Sa to Sc-Im-BCD with the exception of E-SOa, where the 
dust temperature is higher than in star-forming systems proba- 
bly because of the different nature of their dust heating sources, 
c) SPIRE colours can be used to discriminate thermal from syn- 
chrotron emission in radio galaxies, d) In contrast to the warm 
dust, the colour temperature /350//500 index decreases with 
star formation activity and increases with metallicity. This ad- 
mittedly weak evidence might be indicative of an overabundance 
of cold dust or an emissivity parameter/? < 2 in low metallicity, 
active systems. 
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